Ratchet effect for cold atoms in an optical lattice 
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The realization of a directed current for a quantum particle in a flashing asymmetric potential is 
studied. It is found that a positive current, i. e. in the direction expected for a conventional diffusive 
ratchet, can be attained at short times in the limit where the potential is weak and quantum diffusion 
dominates, while current reversal is obtained for stronger potentials. A single parameter, the ratio 
between the kicking frequency and the optical lattice potential strength, determines both the degree 
to which quantum effects dominate, and the possibility of obtaining a ratchet current. The effect 
should be readily observable in experiments. 



The study of nanoscale quantum pumps is of consid- 
erable current interest, motivated by both the quest to 
explore basic quantum phenomena, by recent significant 
experimental progress, and also by the potential practi- 
cal implications in terms of artificial devices Q . Ratchet 
models have been introduced and studied for processes 
that allow extraction of work from diffusive motion with- 
out application of macroscopic gradients. For example, 
ratchets have been discussed in the context of biological 
motors and for manipulating magnetic fiux in artifi- 
cially nanostructured superconductors 

With the advent of optical lattices, the possibility has 
opened up to study dynamics of systems of quantum par- 
ticles in the absence of external noise. These provide a 
novel setting for ratchet effects and interest in purely 
quantum ratchets has therefore grown recently. In rock- 
ing ratchets, where an alternating bias field is imposed 
on the system, current reversal due to quantum effects 
has been observed 0, 0- In an optical lattice, trans- 
port was observed in a setup that utilizes oscillations be- 
tween internal states of the atoms to create an effective 
ratchet potential that alternates between two configura- 
tions ^] . Another realistic experimental setup for atoms 
in an optical lattice is a flashing ratchet, i.e., an asym- 
metric potential that is turned on and off periodically in 
time. Reference carried out a theoretical study of such 
a system in the case of a small effective Planck constant, 
and found that no long-time transport was present in the 
case of strictly periodic flashing. They devised instead a 
scheme with unequal spacing between the flashes, which 
was shown to result in greatly enhanced transport. Ref- 
erence studied a similar system experimentally, but 
without any asymmetry in the potential, and found at 
short times a greatly enhanced spreading in momentum 
space due to quantum resonances at certain parameter 
values (however, because of the geometry, no directed 
transport could be attained). 

In this paper we study a nonlinear quantum pump act- 
ing on a particle trapped in an optical lattice, driven by 
a flashing ratchet mechanism. The focus of the paper 
is to determine the conditions for directed transport in 
a one-dimensional flashing ratchet, in the regime where 
quantum effects are appreciable. The aim is to sort out 



the dependencies on the various parameters of the system 
and to provide physical explanations for the parameter 
dependences and the absence of long-time transport. 

System - Consider a particle in one dimension subject 
to the Schrodinger equation 
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where 



v{x) = sin X + a sin 2x 



(1) 
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is the ratchet potential, which we assume to be composed 
of two standing laser waves. The parameter a controls 
the skewness of the periodic potential that is the origin 
of the ratchet effect. A value of a between and 0.5 
approximates a sawtooth potential where the sawteeth 
lean to the left; this is expected to yield a positive current 
in the usual case of classical diffusive motion. The proflle 
is more asymmetric for larger values of a. The potential 
is assumed to be flashed on and off at periodic intervals. 
To facilitate calculations we have taken the pulses to be 
delta functions in time. The case a = 0, corresponding to 
a pure sine potential, is formally equivalent to the kicked 
rotor and has been extensively studied theoretically and 
experimentally The addition of a second 

sine wave adds complexity to this classic problem and 
opens for the possibility of a ratchet effect j3- 

The units in Eq. are chosen such that the temporal 
period of the flashing and the spatial period of the lattice 
are unity. In terms of physical quantities, the effective 
Planck constant is 
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where T is the period of the flashing and uiji — hk^/2m 
is the recoil frequency of the applied laser fleld, with fci 
the recoil momentum giving a lattice period {2kL)~^ for 
the optical potential felt by the atoms. The effective 
potential strength is 
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where i?oL is the amphtude of the periodic potential of 
the laser beam with the lower frequency. We have implic- 
itly defined the auxiliary quantity P = K/h in Eq. |0J, 
because as we shall see, it determines much of the physics 
of the system. In addition, P has a simple interpretation: 
it is just the ratio between the potential strength and the 
kicking frequency. 

When the potential flashes are assumed to be delta- 
function kicks, the time evolution can in principle be de- 
scribed analytically: in momentum space the wave func- 
tion at time t is given by 
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J2 F{k^q)i,{q,t-l), (5) 
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where the wavenumber difference q — k is restricted to 
integer values because of the periodicity of the potential, 
and the function F is given by 



F{k)= J2 Jk-2p{P)Jp{aP). 



(6) 



p— — oo 



Jn are Bessel functions. The dependency of F on P and a 
is suppressed for brevity. In the kicked-rotor case, a = 0, 
all terms except p ~ vanish in © and there results 
F{k) = Jk{P)- Although the analytical formula, Eqs. 
(j5l6ll . is useful for studying analytical properties, it is 
easier to solve for the time evolution numerically using 
a FFT algorithm. We assume an initially homogeneous 
wave function with zero momentum. The time evolu- 
tion will because of the periodicity bring in only momen- 
tum components V'(^) with integer k. The wave function 
is propagated in time, one period in each step, by first 
multiplying the wave function components in k space by 
a phase factor exp(— then Fourier transforming 
the wave function and applying the phase factor associ- 
ated with the potential, exp(— iP?j(a;)), to the real-space 
wavefunction before Fourier transforming back. 

Effective Planck constant - Figurc^displays the mean 
wave number (fc) of the wavepacket over the first Tav 
temporal periods, as a function of the effective Planck 
constant k. The magnitude of the sin 2a;-term in the po- 
tential is chosen as a = 0.3. The main features of the 
wave number-vs-Planck constant curve do not depend 
critically on the value of a, but as a rule the drift is 
stronger for larger a. We have held P constant, so that 
the potential strength in each of the panels of Fig.^varies 
along the abscissa ds, K = Pk. Consider first the topmost 
panel, where P = 0.5 and Tav = 20. The curve is seen 
to be dominated by a series of resonances, but overall 
there is positive transport in the anticipated direction. 
Shown is only the curve for < S < 47r, because it is pe- 
riodic in k with period 47r. When P is kept constant, the 
potential-energy part of the time evolution operator does 
not depend on k, but only the kinetic term exp(S:A;^/2) 
does, and since spatial periodicity restricts the values of 
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FIG. 1: Mean wave number of the wavepacket over the first 
Tav kicks, as a function of the effective Planck constant fe. 
The asymmetry of the ratchet potential has been chosen to 
a — 0.3. Potential strength P and averaging time Tav are as 
indicated in each panel. 



k to integers, the change k ^ k + An does not affect the 
dynamics. It is thus not obvious that k should be thought 
of as the parameter that determines whether the system 
is classical or quantum mechanical. Moreover, according 
to Eq. (jSJ, changing S to — ^; or 47r — ^ only results in 
an overall complex conjugation of the wave function and 
no change of the physics. This explains the reflection 
symmetry of the curve about the point k = 2tt. 

The resonances in k for integer and half-integer mul- 
tiples of IT are well known from the kicked-rotor case, 
a = H, H E^. The algebra pursued in Refs. [1 E3 
can be carried over to the present, asymmetric case and 
yields similar results. Details of the calculation will be 
presented elsewhere fl^]. 
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Long-time behavior - The second panel in Fig. 
shows the mean wave number over an averaging time 
of Tav = 100 kicks instead of 20. Only the resonances 
remain in this curve; away from these the mean wave 
number has averaged out to nearly zero. Apparently, no 
transport can be expected at long times, except at res- 
onances. This was also found in Ref. and the reason 
is quantum localization. In the kicked rotor problem, it 
is well known that a quantum system will at short times 
undergo diffusion in momentum space, just like in the 
classical case, at a diffusion rate D = K'^/2. However, 
after a so-called break time t* oc D/k^, the quantum time 
evolution ceases to follow the corresponding classical evo- 
lution and the momentum distribution does not broaden 
any more. This is analogous to Anderson localization, al- 
though the phenomenon takes place in momentum space 
and not in real space The localization persists also 
in the case of an asymmetric potential, a ^ 0, and seems 
to forbid a persistent current (see, however, Ref. 0). Ob- 
serve that with our definitions, the break time is simply 
t* oc P^. In other words, a large P means a more pro- 
nounced quantum-classical correspondence. The relevant 
parameter determining the quantum/classical character 
of the dynamics is thus P and, perhaps somewhat coun- 
terintuitively, not k. 

The present parameter value, P = 0.5, is far into 
the quantum regime since it results in a break time 
t* less than one, and thus no quantum-classical cor- 
respondence can be expected. Figure HI displays the 
time evolution of the wave packet for a representative 
choice of parameters away from resonances: k — 1.7, 
and as before, P — 0.5 and a — 0.3. The lower panel 
of Fig. 121 shows the momentum asymmetry, defined as 
the deviation of the density from its symmetrized value: 
6n{k) = n{k) — {n{k) + n{—k))/2. It is seen that the den- 
sity distribution initially broadens, with a slight tendency 
towards positive fc, whereafter the width of the distribu- 
tion saturates due to localization, because of the short 
break time. From that point, the momentum distribu- 
tion begins to oscillate back and forth within the wave 
packet, yielding a zero net momentum at longer times. 
The initial tendency towards positive k is in accordance 
with our expectancies for a ratchet; this is further dis- 
cussed below. 

Kick strength - As discussed above, it is customary 
to discuss the physics of the present system in terms of 
the kick strength K, because the Schrodinger equation 
(P) assumes the intuitively expected form, and more- 
over, the effective Planck constant k assumes the role 
of a quantum parameter. However, as we have seen, 
P — K/k is in fact the natural parameter that deter- 
mines the degree to which the time evolution is quantum 
mechanical, and furthermore, as we now show, it turns 
out to determine whether a short-time ratchet effect is 
present. In the third panel of Fig. ^ we report the result 
for slightly stronger kicks than in the examples above. 
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FIG. 2: Time evolution of the wave function in momentum 
space. Upper panel displays the density distribution; the 
logarithm is taken in order to enhance visibility. Each ver- 
tical cross-section of the plot represents the squared wave 
function n{k,t) = \Tp{k,t)f at a time instant t. The lower 
panel shows the asymmetry in the momentum distribution, 
Sn{k,t) = n{k,t) - {n{k,t) +n{-k,t))/2. 
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FIG. 3: Regimes of positive (white areas) and negative trans- 
port (black areas) for the quantum ratchet with a = 0.3. The 
mean momentum is taken over the first Tav = 20 periods. The 
first and third panels in Fig. correspond to the horizontal 
cross-sections of the present curve at constant P — 0.5 and 
P — 2, respectively, while the middle panel of Fig. 0is taken 
at a different Tav. 



P — 2. The average is taken over twenty periods. For 
this kick strength, the curve has a considerably more bro- 
ken appearance. In particular, the drift is no longer pos- 
itive for all values of k: it is just as likely to be negative. 
Indeed, we do find that for small values P ^ 1, the drift 
is always positive, except at anti-resonances, but it can 
be negative for larger P. Figure |21 shows the regions in 
P-k space where the drift over the first twenty kicks is 
positive and negative, respectively. It is seen that for 
large values of P or X and for < ^ < tt, the black 
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regions run along lines of constant K, indicating that in 
the classical regime, the dynamics depends on K rather 
than P. On the other hand, what is interesting here is 
the dividing line below which the drift is purely positive, 
and that is clearly a line of constant P, not K. 

The observation of negative transport (current rever- 
sal) for quantum ratchets in the zero-temperature limit 
has in the case of a rocked ratchet been explained with a 
tunneling argument 0. That argument is not applicable 
for a delta- function flashing ratchet; obviously another 
mechanism is at play here. Consider the time evolution 
of a wavepacket originally concentrated within a well. 
Between two kicks, it undergoes diffusion according to 
the Schrodinger equation (QJ and diffuses out to a final 
position uncertainty AaJdiff ~ S in the present units. The 
particle also receives an impulse from the periodic poten- 
tial, which is of the order Ap ^ K \i the particle initially 
resided on one of the slopes of the potential, and hence 
travels a distance comparable to AaJbaiUstic ~ if in one 
period. We now see that the condition for a ratchet effect 
at short times, 

P<1, (7) 
translates into a condition on the physical time scales, 

Axballistic ^ Axdiff. (8) 

In summary, we have investigated the conditions for 
directed transport for a quantum particle in a flashing 
ratchet potential. It is found that at short times, a di- 
rected current will flow in the anticipated direction if the 
kick strength P, defined as the ratio between the ampli- 
tude of the potential and the flashing frequency, is less 



than unity; this coincides with the limit in which the 
motion is inherently quantum mechanical, and quantum 
diffusion dominates over ballistic transport. It would be 
interesting to look for these findings in experiment: the 
levels of precision, noise and parameter control attain- 
able in current experiments 0, 0] are certainly sufficient 
to test the predictions made in the present paper. 
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